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ABSTRACT 


A  method  of  calibrating  an  electromagnetic 
seismometer  using  a  variation  in  the  surrounding 
magnetic  field  has  been  developed.  A  pair  of  coils, 
mounted  externally,  connected  in  a  series-opposing 
configuration  are  used  to  apply  a  step  in  force  to 
the  magnet  suspended  within  a  Willmore  Mk  II  seismo¬ 
meter.  The  force  is  determined  from  a  measurement  of 
the  steady  state  deflection  of  the  magnet  with  a  linear 
variable  differential  transformer.  The  output  signal 
is  digitized  and  recorded  on  a  nine-track,  800  bpi  IBM 
compatible  magnetic  tape.  A  z-transform  is  used  to 
transform  the  digitized  data  into  the  frequency  domain. 
This  is  then  differentiated  and  normalized  so  as  to 
give  the  velocity  sensitivity  of  the  seismometer. 

This  method  provides  a  rapid  and  accurate  method 
of  calibrating  the  seismometer.  Once  the  coils  are 
mounted  no  contact  is  needed  with  the  seismometer  so 
that  calibrations  can  be  done  regularly  and  frequently 
even  on  seismometers  which  are  not  easily  accessible. 
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CHAPTER  1 


CALIBRATION  PROCEDURES 

1 . 1  Introduction 

Seismic  studies  involve  more  than  the  arrival 
times  and  directions  of  the  prominent  waves.  They 
require  a  knowledge  of  the  amplitude  and  frequency 
content  of  the  ground  motion.  It  is  necessary,  then, 
to  determine  accurately  the  characteristics  of  the 
system.  This  requires  a  calibration  technique  to  deter¬ 
mine  the  magnification  or  velocity  sensitivity  and  the 
phase  response  of  the  seismometer. 

There  are  several  methods  employed  in  obtain¬ 
ing  the  frequency  response  of  the  seismometer  from  its 
output.  Espinosa  et  al  (1962)  used  a  transient  technique 
which  involved  taking  the  ratio  of  the  Fourier  trans¬ 
forms  of  the  output  and  the  input  transients  to  yield 
both  the  phase  and  amplitude  response  of  the  system  in 
the  frequency  domain.  The  input  transient  used  was 
either  a  step  function  or  a  delta  function.  Mitchell 
and  Landisman  (1968)  showed  that  the  direct  Fourier 
analysis  gives  poor  results  when  microseisms  are  pre¬ 
sent  and  are  not  negligible  in  size  compared  to  the 
output  resulting  from  the  calibration  transient.  They 
suggested  that  a  least  squares  fit  be  used  on  the  seis¬ 
mometer  output  in  the  time  domain  to  determine  the 
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parameters  of  the  seismometer.  The  frequency  response 
can  then  be  calculated  using  the  transfer  function 
mathematically  determined  for  the  system.  Certain 
calibration  techniques  enable  the  frequency  response 
to  be  determined  directly  from  the  output  by  using  a 
sinusoidal  input  at  the  desired  frequencies  (Murphey 
et  al  1952).  Though  less  is  involved  in  processing  data 
in  this  form,  this  method  is  more  time  consuming  than 
the  others  because  it  requires  readings  at  several  fre¬ 
quencies  for  accurate  analysis. 

The  advent  of  the  computer,  improved  seismometer 
design  and  the  use  of  digital  recording  methods  (Burke 
et  al  1970)  has  given  the  seismologist  more  accurate 
methods  of  handling  data.  It  has  introduced  the  use  of 
an  active  seismograph  system  (use  of  an  amplifier  and 
digital  recorder  instead  of  a  galvanometer)  which  permits 
the  seismologist  to  alter  the  response  of  the  seismometer 
system  with  feedback  controlled  amplifiers  and  filters  to 
get  the  desired  response  from  the  system.  This  has 
allowed  seismologists  to  attempt  studies  of  recorded 
pulse  shapes  as  modified  along  the  propagation  path,  a 
process  which  requires  an  accurate  knowledge  of  the  seis¬ 
mometer  transfer  function. 

Deviations  of  the  system  parameters  from  assumed 
values  can  lead  to  serious  errors  in  phase  and  amplitude 
corrections.  The  fact  that  these  parameters  may  vary  due 
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to  various  influences  such  as  temperature  and  pressure 
or  fatique  of  the  material  of  the  structural  element, 
and  that  these  variations  may  occur  aperiodically  or 
systematically,  over  a  long  period  of  time  (one  year) 
or  a  short  period  of  time  (one  day)  necessitates  regular 
and  frequent  calibration  of  the  seismometer. 

1 . 2  Calibration  Techniques 

There  are  basically  three  types  of  seismometers: 
the  moving-coil  or  moving-mass  type,  the  changing  flux 
type  and  the  magneto-striction  type.  Shima  (1960) 
pointed  out  that,  in  earthquake  analysis,  the  third  type 
is  not  used  and  the  second  type  is  used  to  observe  only 
minute  and  distant  earthquakes  because  of  their  non-linear 
response  to  strong  earth  motions.  Most  seismometers  are 
of  the  first  type.  It  is  necessary  that  a  calibration 
technique  be  applicable  on  this  type  of  instrument  on  a 
routine  and  frequent  basis. 

a)  Weight  Lifting 

Shima  (1960)  proposed  a  simple  method  for  deter¬ 
mining  the  voltage  sensitivity  of  a  seismometer  (figure 
1-1).  By  passing  a  current  (I)  through  the  coil,  a 
force  (AI ,  where  A  is  the  motor  factor  of  the  seismo¬ 
meter)  is  exerted  on  the  seismometer  mass  (M)  and  dis¬ 
places  it.  The  seismometer  is  drawn  back  up  to  its 
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Figure  1-1  Determining  the  voltage  sensitivity  of  a 

seismometer  by  weight  lifting  technique 
(Shima,  1960) . 


original  position  by  adding  a  mass  (m)  to  the  left  hand 
side  of  the  pan  balance.  A  string  attached  to  the  right 
hand  side  of  the  pan  balance  then  pulls  the  seismometer 
up  to  its  original  position  which  is  determined  opti¬ 
cally.  The  value  of  A  can  then  be  determined  from  the 
conservation  of  force  which  requires  that  mg  =  AI  where 
g  is  the  acceleration  of  gravity. 

This  method  provides  only  one  of  the  instrument 
constants.  A,  and  so  the  frequency  response  of  the 
system  cannot  be  determined.  It  is  also  impractical 
for  frequent  field  calibrations  as  it  requires  access 
to  the  seismometer . 

b)  Shake  Table 

In  order  to  obtain  the  frequency  response  of  the 
system  it  is  necessary  to  apply  a  transient  and  observe 
the  resulting  oscillations  of  the  seismometer  mass  or 
apply  a  sine-wave  and  record  the  output  amplitude. 

One  way  of  doing  this  is  with  a  shake  table.  This  is 
a  vibrating  platform  which  can  be  displaced  a  given 
amount  at  a  given  frequency  and  have  a  constant  dis¬ 
placement  over  a  wide  range  of  frequencies.  Accurate 
shake  table  calibrations  require  an  accurate  knowledge 
of  the  frequency  and  magnitude  of  the  shake  table  vibra¬ 
tions.  Improved  shake  table  accuracy  has  been  brought 
about  by  the  use  of  piezoelectric  crystals  (Hoover  et  al 
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1965;  T.  Bardeen  1950)  as  a  driving  mechanism  by  which 
voltage  input  is  transformed  to  motion.  The  transient 
input  of  the  shake  table  can  be  an  impulse,  step  or 
sinusoidal  wave  (calculations  -  Hoover  et  al  1965) . 

Shake  table  calibrations  cannot  be  easily  used 
for  frequent  field  calibrations  in  that  they  demand 
contact  with  the  seismometer  which  may  be  buried  to 
reduce  temperature  and  pressure  variations.  The  shake 
table  is  also  too  bulky  to  be  used  in  field  calibrations. 
It  is  a  useful  process  by  which  to  initially  calibrate 
the  seismometer  but  it  needs  to  be  paired  up  with 
another  method  of  calibration  for  field  purposes. 
Difficulty  also  arises  in  using  the  shake  table  for 
the  calibration  of  horizontal  seismometers  in  that  a 
slight  tilt  in  the  horizontally  vibrating  shake  table 
will  result  in  the  addition  of  an  external  force  on  the 
system  due  to  gravity. 

c)  Weight  Adding  or  Removing 

Another  technique  of  seismometer  calibration  is 
to  add  a  small  mass  to  or  remove  it  from  the  moving  mass 
of  the  seismometer.  The  removal  or  addition  of  a  smaller 
mass  m  is  the  equivalent  to  having  a  force  mg  acting  on 
the  seismometer.  If  the  removal  or  addition  is  done 
instantaneously  it  is  equivalent  to  displacing  the  mass 
a  distance  Ax  =  mg/k  (where  k  is  the  spring  constant  of 
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the  seismometer)  from  its  rest  position  and  releasing 
it.  As  this  procedure  requires  contact  with  the  moving 
mass  of  the  seismometer  it  is  not  useful  for  field 
purposes . 

A  variation  on  the  removal  of  the  mass  was 
suggested  by  Dix  (1952) .  He  proposed  the  seismometer 
be  placed  on  a  platform  supported  by  a  spring  having  a 
considerably  longer  natural  period  of  oscillation  than 
the  seismometer  and  that  the  platform  have  a  mass  sus¬ 
pended  from  it  (figure  l-2a) .  A  sudden  removal  of  the 
suspended  mass  results  in  the  seismometer  experiencing 
a  force,  mg,  which  can  be  considered  a  step  where  the 
period  of  the  supporting  spring  is  considerably  larger 
than  the  length  of  the  record  (figure  l-2b) .  The 
necessity  of  the  suspended  platform  makes  this  a  poor 
procedure  for  field  calibrations. 

d)  Maxwell  Bridge 

The  seismometer  is  an  electromechanical  device. 
Most  procedures  thus  far  have  used  the  mechanical  aspect 
of  the  seismometer  for  calibration  purposes.  Another 
type  of  procedure  is  to  apply  a  force  to  the  moving 
mass  of  the  seismometer  by  putting  a  current  through 
the  seismometer  coil.  Putting  a  d.c.  current  through 
a  seismometer  coil  to  displace  the  mass  and  then  by  means 
of  a  switch  simultaneously  shutting  off  the  current  and 
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Figure  l-2a  Apparatus  used  in  the  calibration  of  the 

seismometer  using  step  in  acceleration 
(Dix,  1952) . 


Figure  l-2b 


Force  acting  on  seismometer  due  to  the 
removal  of  mass  (m)  as  a  function  of 
time  (Dix  ,  1952 )  . 


connecting  the  seismometer  to  a  recording  device 
(current  release  test:  K.G.  Barr  1964;  Hoover  et  al 
1965)  has  been  used  as  a  method  of  field  calibration 
to  complement  other  methods  of  seismometer  calibration. 

Willmore  (1959)  introduced  the  use  of  a  Maxwell 
bridge  to  input  the  current  into  the  seismometer 
(figure  l-3a) .  His  calibration  procedure  consisted  of 
four  steps: 

1.  Balancing  the  bridge  with  the  seismometer 
clamped. 

2.  Reading  the  size  of  the  deflection  of  the 
galvanometer  at  several  frequencies  with  the 
seismometer  unclamped. 

3.  Amplitude  readings  are  taken  at  several  fre¬ 
quencies  with  the  seismometer  clamped  and  the 
bridge  on  substitution  input. 

4.  Repeat  step  3  with  the  seismometer  unclamped. 

Any  current  put  through  the  coil  exerts  a  force, 

Ai  ,  (A  -  motor  factor,  i  -  input  current)  on 
c  c 

the  seismometer  mass.  The  signal  on  the  measuring  device 
(a  galvanometer  in  Willmore' s  paper)  is  then  greater  than 
that  due  to  the  input  current  by  an  amount  due  to  the 
resulting  oscillations  of  the  seismometer  mass. 

Balancing  the  bridge  reduces  the  input  current,  as  seen 
by  the  measuring  device,  to  zero.  The  observed  motion 
is  then  entirely  that  corresponding  to  a  force,  Aic, 
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Figure  l-3a  Maxwell  bridge  used  in  seismometer 

calibration  (Willmore,  1959). 


Figure  l-3b  Modified  Maxwell  bridge  used  in 

seismometer  calibration  (Shima  et  al , 
1966)  . 


acting  on  the  seismometer  mass.  This  is  equivalent  to 

a  ground  acceleration  of  -Ai  /M  (M  is  the  seismometer 

c 

mass) .  The  substitution  input  enables  an  absolute 
calibration  to  be  made.  Though  Willmore' s  original 
measuring  device  was  a  galvanometer,  an  oscilloscope 
gives  a  simpler,  more  accurate  and  direct  method  of 
calibration  (White  1970). 

A  transient  input  in  the  form  of  a  step  function 
as  opposed  to  the  sine-wave  used  by  Willmore  -  was  used 
by  Espinosa  et  al  (1962)  to  determine  the  frequency 
response  of  a  long  period  seismometer  by  comparing  the 
Fourier  transforms  of  the  output  and  input  signals.  At 
the  low  frequencies  considered,  the  inductive  reactance 
of  the  coil  can  be  ignored  and  the  seismometer  merely 
regarded  as  a  resistor  which  has  no  frequency  dependence 
The  Maxwell  bridge  can  then  be  replaced  by  a  Wheatstone 
bridge.  However  at  higher  frequencies,  1  to  20  cps,  the 
seismometer  cannot  be  represented  by  just  a  resistor  or 
as  a  resistor  in  series  with  an  inductance  as  done  by 
Willmore.  No  single  balance  point  can  be  found  unless 
losses  caused  by  hysteresis  and  eddy  currents  are  consi¬ 
dered.  Shima  et  al  (1966)  took  this  into  account  by 
representing  the  seismometer  by  a  resistance  in  parallel 
with  an  inductance  and  both  in  series  with  a  second 
resistance.  Their  modified  Maxwell  bridge  (figure  l-3b) 
can  then  be  balanced  at  one  point  for  frequencies  up  to 


- 
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20  cps  and  the  calibration  carried  out  using  a  step 
input . 

The  present  calibration  procedure  at  the  Depart¬ 
ment  of  Physics,  University  of  Alberta  takes  into 
account  hysteresis  and  eddy  currents  by  using  sine- 
waves  of  various  frequencies  as  input  and  balancing 
the  bridge  at  every  frequency.  It  has  been  found  that 
this  technique  is  not  suitable  for  field  calibrations 
in  that  it  requires  contact  with  the  seismometer  for 
clamping  and  unchamping  in  addition  to  being  a  long  and 
tedious  method. 

e)  Calibration  Coil 

The  electric  current  can  also  be  applied  to  a 
separate  calibrating  coil  mounted  in,  on  or  around  the 
seismometer.  Stewart  et  al  (1967)  proposed  an  external 
coil  for  the  calibration  of  a  Willmore  Mk  II  seismometer 
which  was  mounted  on  top  of  the  seismometer.  The 
positioning  rod  of  the  seismometer  is  replaced  by  a 
connecting  rod  with  a  small  magnet  mounted  on  the  end 
of  it;  mass  is  kept  constant.  A  current  passing  through 
the  calibrating  coil  exerts  a  force  on  the  seismometer 
mass.  The  type  of  forcing  function  acting  on  the  seis¬ 
mometer  depends  on  the  signal  put  through  the  coil.  The 
motor  function  of  the  calibration  coil  is  determined  by 
the  weight-lifting  method. 
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Certain  other  seismometers  have  an  auxiliary 
coil  built  into  them  for  calibration  purposes.  How¬ 
ever  this  is  not  always  the  case.  The  method  of 
calibration  in  this  thesis  is  the  use  of  two  coils 
mounted  externally  to  the  seismometer  with  the  current 
through  one  coil  flowing  in  the  opposite  direction  to 
the  current  in  the  other  coil.  The  two  coils  are 
positioned  so  as  to  minimize  the  induced  electrical 
signal.  A  step  function  in  current  through  the  coils  which 
changes  the  magnetic  field  in  the  volume  of  the  seismo¬ 
meter  acting  on  the  suspended  magnet  is  equivalent  to 
a  step  in  the  inertial  field  acting  on  the  seismometer 
mass.  The  motor  function  of  this  arrangement  is  deter¬ 
mined  by  measuring  the  resulting  displacement  of  the  mass, 
A£,  using  a  linear  variable  differential  transformer 
mounted  on  the  end  of  the  positioning  bar  on  the  Willmore 
Mk  II  seismometer.  The  motor  function  could  also  be 
determined  by  initially  calibrating  the  seismometer  using 
another  technique.  This  particular  method  is  feasible 
for  frequent  field  calibrations  in  that  it  does  not  re¬ 
quire  contact  with  the  seismometer  and  can  quickly  and 
easily  be  carried  out.  A  sine-wave,  step  function  or 
impulsive  function  may  be  put  through  the  coil. 


;irc 


CHAPTER  2 


THE  SEISMOMETER  TRANSFER  FUNCTION 
2 . 1  Transfer  Function 

The  output  of  the  seismometer  to  an  input  of  a 
delta  function  in  ground  displacement  or  velocity  is 
known  as  the  displacement  or  velocity  transfer  function 
respectively.  The  expression  for  the  transfer  function 
can  be  derived  either  from  an  electrical  equivalent  net¬ 
work  of  the  complete  seismograph  system  (Kollar  et  al, 
1966)  or  a  system  of  differential  equations  describing 
the  elements  of  the  seismograph  system,  as  presented  in 
this  chapter. 

The  seismometer  used  was  a  Willmore  Mk  II  seismo¬ 
meter  (Willmore  et  al ,  1963).  It  consists  of  a  magnetic 
mass  suspended  over  a  coil  by  a  spring  system.  This  is 
represented  in  figure  2-1  with  damping  represented  by  the 
dash  pot.  x  is  the  displacement  of  the  seismometer  case, 
E,  the  displacement  of  the  suspended  mass. 

When  the  frame  is  in  motion  upwards  (positive  x) 
the  forces  on  the  mass  M  are 

1)  an  elastic  force  (k  times  the  net  extension  of  the 
spring;  k  is  the  elastic  constant;  &  is  the  rest  exten¬ 
sion  of  spring  due  to  seismometer  mass) 

F  =  kU  +  x  -  £)  ; 

2)  the  damping  force  (retarding  a  times  the  velocity 
of  the  mass  with  respect  to  the  frame;  a, the  damping 
constant  is  the  sum  of  viscous  and  electric  damping) 
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Figure  2-1  The  coordinate  system  of  the  seismometer. 

Arrows  indicate  positive  displacement. 
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F  =  -a  (  £  -  x) 
a 

3)  force  of  gravity  (acting  downwards) 

F  =  -Mg 

4)  magnetic  force  due  to  current  through  external 
coil 

m 

The  sum  of  these  forces  causes  an  acceleration 

•  • 

of  the  mass  £.  Newton’s  Second  Law  then  requires 

k(&  +  x  -  £)-a(£-x)-Mg  +  F^  =  M? 

Let  y  =  x  -  £  be  the  relative  displacement  of  the 
frame  and  mass 


k£  +  ky  +  ay  -  Mg  =  -F  +  M(x-y)  . 


k£  is  the  original  extension  of  the  spring  due  only  to 
the  force  of  gravity  on  the  mass  M.  Therefore 

k£  =  Mg 


My  +  ay  +  ky  =  Mx  -  F 


m 


(2.1) 


2  k 

a)  *=  —  U  =  natural  period  of  the  seismometer  in 
n  M  n  c 


2X  =  a 

a)  k 

n 


radians ) 


(£  =  damping  factor) 


y  +  — ^  y  +  -X?  =  —77  -  -7^-  (equation  of  motion) 

“n  u2  co2  k 
n  n 


' 
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Case  Is  No  magnetic  force  (i.e.  Fm  =  0) 


Taking  the  Laplace  transform  of  the  equation  of 


motion 


2  2 

Y(s)  +  •  s  •  Y(s)  +  ^  Y(s)  =  X(s)  . 

n  co  go 

n  n 

The  seismometer  is  a  physically  realizable  instrument  so 
the  Laplace  transform  can  be  used  to  solve  the  differ¬ 
ential  equation.  The  Laplace  transform  of  y(t)  is  Y(s). 


Normalize  this  expression  using  S  = 


co 


n 


Y(S)  = 


1  +  2  CS  +  S' 


X(S) 


(2.2) 


the  transfer  function  is  then 


Y  (s ) 
X  (s ) 


1  +  2  CS  +  S 


2  =  VS) 


for  a  delta  function  in  ground  displacement  X(S)  =  1 
and  the  transfer  function  is  represented  by  Yq(S). 

The  output  voltage  of  the  seismometer  is  repre¬ 
sented  by 

e(t)  =  -Ay(t)  (A  =  motor  factor) 


taking  the  transform 


E  (s )  =  -As  Y  (s ) 

E(S)  =  -A  cor  S  Y  (  S)  .  (2.3) 

If  x(t)  is  the  input  displacement  then  x(t)  is 
the  input  velocity 
v  (t)  =  x  (t) 


V  (s )  =  s  X  (s ) 


' 
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or  V  (S)  =  a)  S  X  (S) 

n 


(2.4) 


substituting  equations  (2.3) 

(2.2) 

E(S)  _  S 

-Ago  (S) 


n 


1  +  2?S  +  S' 


and  (2.4) 
V(S) 

03nS 


into  equation 


E  (S) 
V(S) 


-AS' 


1+  2?S  +  S 


2  ~  Eo(S) 


(2.5) 


Eq(S)  is  the  velocity  transfer  function  -  the  transform 
of  the  output  voltage  for  a  delta  function  input  in 
velocity . 


Case  2:  No  ground  motion  and  a  general  function  for  the 
magnetic  field  motion 

The  magnetic  force  Fm  can  be  represented  by 

=  kD(t)  where  k  is  the  spring  constant  and  D(t)  is 

the  displacement  of  the  mass  from  its  rest  position  due 

to  the  force  F  . 

m 

The  equation  of  motion  is  then 


y  +  7T  y  +  ^2  =  -D(t) 

n  ^ 


or  transformed  and  normalized 


Y (S)  [1 +  2£S  +  s  ]  =  -D  (S) 


-D(S) 


1  +  2  CS  +  S‘ 


Y(S) 


- 

' 
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expressing  Y(S)  in  terms  of  voltage 


E  (S)  =  Aoj  S  Y  (S ) 
n 


therefore 


D 

o 


(S) 


E  (S) 
D(S) 


+A00  S 

_ n _ 

1  +  2CS  +  S2 


(2.6) 


where  Dq(S)  is  the  transfer  function  due  to  the  magnetic 
force.  The  velocity  transfer  function  is 


Eq  (S)  = 


■AS 


1  +  2  £S  +  S' 


Therefore  the  velocity  transfer  function  can  be  derived 
from  the  magnetic  transfer  function. 


E  (S)  =  -  —  D  (S)  . 

O  h)  O 


2 . 2  A  Step  in  Input 

Take  the  case  where  the  magnetic  force  is  due 
to  a  step  in  the  magnetic  field  affecting  the  seismometer. 


F  =  kD (t )  t  >  0 

m 

F  =  kA&  t  £  0 

m 

=  0  t  <  0 

.*.  D  (t )  =  A£  t  >  0 

=  0  t  <  0 


A  £ 


D  (S) 


AZ 

s 


a)  S 

n 


•- 
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Then  we  get  from  (2.6) 


mag 


(S)  = 


+AA& 


1  +  2  £S  +  S' 


(2.7) 


the  velocity  transfer  function  is  then 


Eo(S)  = 


~  Em  (S) 
/a  36  rQ8.g 


The  velocity  transfer  function  can  then  be  obtained 
from  a  step  in  the  magnetic  field  by  multiplying  the 


transfer  function  of  the  output  data  by 


(o. 


S  =  -  j 


0) 


n 


A£ 


S2  = 


60 


n 


E0(s)  = 


+60 


60. 


n 


j-r-  E  (S) 

A  ic  mag 


The  value  for  A£  can  be  obtained  by  measurement  of 
the  output  of  the  linear  variable  differential  transformer 


2 . 3  A  Sine-wave  Input 

The  transfer  function  can  be  obtained  directly, 
by  using  a  sine-wave  input  in  the  magnetic  field.  In 
this  case  the  equation  of  motion  can  be  rewritten 


y  (t)  +  —  y  (t)  +  -ij-  y  (t)  =  -A£  sin  oo1  (t)  .  (2.8) 

wn  60  ^ 

n 


as 
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The  solution  of  this  differential  equation  is 
of  the  form 


y(t)  =  yp(t)  +  yh(t) 


y  (t)  =  k,  sin  03,  (t)  +  k~  cos  oin  (t) 

J  p  1  1  2  1 


substituting 

back  into 

(2.8)  to  solve 

for 

i — 1 

44 

and  k2 

03 

2C  — 

03 

i — 1 

44 

sin  03^ 

2 

03 

(t)  -  (1-  —■) 

k2 

sin 

“i 

(t) 

n 

wn 

^1 

+  2C  — 

03 

n 

k2 

COS  03  ^ 

°°1 

(t)  +  (i  -  ~4) 

03Z 

i — 1 

44 

cos 

(t)  = 

n 


Solving  for  k^  and  k2 


y  (t)  = 

P 


“l 

2C  —  Ail 

03 

n 


-Ail  sin  o)^  (t ) 


03 


i  2  “i  2 

<1+^>  +<2?— > 

03  n 

n 


cos  03  ^  (t)  + 


03 


-Ail  (1 - x) 

< 


03^  2  0)  1 

(1-  -4)  +  (2c  — — ) 

03 2  wn 

n 


—  sin  0)  (t )  = 

2  1 


Ail 


2  “i  2 

d-4>  +  (2S  sr> 

03  n 

n 


sin  (o)^t 


-  6) 


where 
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03- 


0  =  arctan 


03 


n 


03 


1  - 


03  ' 

n 


-03  £t 

yh  (t)  =  k3e  n  sin  03r 


-oj  tt 


A-  -2  ■  ■  -  "n 


£  t  +  k^e 


COS  0) 


>/^7t 


but  as  we  are  interested  in  only  the  steady  state  solu¬ 
tion  and  in  that  case  y^(t)  =  0 


y  (t)  =  yp(t)  = 


■AZ 


2  u.  2 
(1  -  4>  +  <2?  ^r> 

03  n 

n 


sin  (o3^t 


-  0) 


The  output  from  the  seismometer  is  then 


esin<“l)  =Ay(t)  ” 


-A£o3^A 


COS  ( 03 2 1 


-  0) 


We  are  interested  in  the  velocity  transfer  function  (2 . 5 ) 
which  is 

,2 


E  ( s )  = 


AS 


1  +  2CS  +  S' 


using  the  substitution 


0). 


S  =  -  j 


03 


n 


■ 
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(0 


CO 


Eo(wl} 


1 

2 

n 


1  - 


W1 

2  j  —  c 

wn 


2 

^1 

2 


co 


n 


solving  for  the  amplitude  and  phase  of  above 

2 


co- 


Amp 


-/ 


2  2 
Real  +  Imag  = 


co 


n 


real 


co 

1 - — 

2 

% 


The  peak  to  peak  output  of  the  sine-wave  input 


is 


e  .  ( co )  = 

sm  1 


-2A£co1A 


and  the  phase  delay  of  the  ouput  from  the  input  is 


cos  (co^t  -  0)  =  -sin(co^t  -  0  -  90) 


Phase  =90+0 


The  velocity  transfer  function  can  be  obtained 

by  multiplying  the  output  amplitude  of  each  value  of 

2 

the  sine-wave  by  co^/2AJlcon  and  the  phase  by  subtracting 
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90°  from  the  measured  phase  difference. 

The  output  of  the  linear  variable  differential 
transformer  in  the  sine-wave  case  is 

V (t)  =  Ly ( t ) 


V  (t )  = 


-  A£L 


2  2 

% 

2 ' 

co 

n 


(Of 

+  (2S 


sin  (co^t 


-  0) 


The  peak  to  peak  value  is 


V(t)  = 


+2A&L 


which  approaches  the  value  2A£L  as  co^->  0  and  since  L 
is  known  A£  can  be  obtained  and  the  velocity  transfer 
function  calculated. 


CHAPTER  3 


THE  CALIBRATING  PROCEDURE 

3 . 1  Introduction 

The  calibrating  procedure  is  to  apply  a  known 
magnetic  force  to  the  suspended  magnet  (mass)  of  the 
seismometer.  The  suspended  magnet  (figure  3-1)  consists 
of  an  open-ended  cylinder  of  Magloy  which  is  a  permanent 
magnet,  a  magnet  keep  which  makes  up  the  inner  pole  piece, 
and  an  outer  pole  piece.  This  geometry  results  in  a 
radial  magnetic  field  across  the  air  gap  in  which  the 

seismometer  coil  moves.  The  relative  motion  of  the 
suspended  magnet  to  the  coil  which  is  fixed  to  the  seis¬ 
mometer  case  induces  an  e.m.f.  in  the  coil.  This  induced 
e.m.f.  is  the  seismometer  output  signal. 

3 . 2  The  Calibrating  Coils 

An  external  magnetic  field  is  applied  to  the  sus¬ 
pended  magnet  by  putting  a  current  through  a  pair  of 
coils  mounted  outside  the  seismometer.  For  experimental 
purposes  the  coils  are  supported  on  a  frame  about  the 
seismometer  (figure  3-2)  though  in  normal  operation  the 
coils  will  be  assembled  into  grooves  machined  in  the 
seismometer  case.  This  will  make  the  coils  an  integral 
part  of  the  case  and  eliminate  the  possibility  of  error 
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magnet 


outer  pole 
piece 


seismometer 

coil 


ENIA  MILD  STEEL 


MAGLOY 


mass  =  4.75  kg 

note:  the  magnet  is  cylindrically  symmetric 


Figure  3-1  The  suspended  mass  (magnet)  in  a  Willmore 

Mk  II  seismometer. 
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Seismometer 


Calibration 

Coils 


Figure  3-2  Calibration  coils  mounted  on  a  frame  about 

the  Willmore  Mk  II  seismometer. 
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due  to  the  dependence  of  input  force  on  the  vertical 
position  of  the  coils  with  respect  to  the  seismometer 
mass . 

Each  coil  consists  of  one  hundred  turns  of  number 
twenty-six  A.W.G.  copper  wire  and  has  a  resistance  of 
seven  ohms.  The  coils  are  wound  on  bobbins  which  have 
an  inside  diameter  of  seventeen  and  a  half  centimeters 
and  depth  of  six  millimeters.  The  two  coils  are  connec¬ 
ted  in  series  and  the  calibration  signal  is  applied  from 
a  Hewlett-Packard  Model  20 3A  function  generator.  A  one 
kilohm  resistor  is  connected  in  series  with  the  generator 
so  that  it  appears  as  a  current  source  to  the  load.  This 
reduces  the  effects  of  the  calibration  coil  parameters 
(inductance)  on  the  input. 

The  two  coils  were  originally  set  up  in  a  Helmholtz 
configuration.  The  coil  separation  was  half  the  inside 
diameter  of  the  coils  and  the  coils  were  connected  in  a 
series  aiding  configuration.  A  change  in  magnitude  of 
a  near  uniform  magnetic  field  was  then  used  to  apply  the 
change  in  force  on  the  suspended  mass.  The  resulting  de¬ 
flections  were  small  resulting  in  output  voltages  of  less 
than  one  volt  per  ampere  of  current  applied  to  the  coils  and 
large  induced  e.m.f.'s  were  observed  in  both  the  seismometer 
coil  and  the  linear  variable  differential  transformer  (LVDT ) . 

The  direction  of  current  flow  in  one  coil  was  then 
reversed  so  that  the  coils  were  now  connected  in  a  series 


- 
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opposing  configuration.  The  deflection  of  the  seismometer 
in  this  configuration  was  greater  by  several  orders  of 
magnitude.  The  mechanism  involved  in  producing  this 
effect  can  be  seen  by  examining  the  magnetic  field 
between  the  two  coils  and  considering  that  the  transla¬ 
tional  force  acting  on  a  magnet  in  a  magnetic  field  depends 
directly  on  the  gradient  of  the  magnetic  field  in  the 
direction  of  the  force  and  the  component  of  the  magnetic 
dipole  (which  represents  the  suspended  magnet)  that  lies 
parallel  to  this  gradient.  The  field  in  the  series  aiding 
configuration  is  a  fairly  uniform  field  and  so  very  little 
net  force  acts  on  the  magnet.  The  series  opposing  confi¬ 
guration  gives  a  large  gradient  in  the  axial  direction 
(illustrated  qualitatively  in  figure  3-3) .  A  larger  force 
then  acts  on  the  magnet.  Changing  the  current  direction 
in  both  coils  changes  the  sign  of  the  gradient  and  the 
force  acts  in  the  opposite  direction. 

The  non-uniformity  of  the  direction  of  the  magnetic 
field  (it  is  not  all  in  the  axial  direction)  enables  a 
force  to  act  in  the  axial  direction  even  if  the  magnetic 
moment  of  the  suspended  magnetic  mass  lies  in  a  radial 
plane.  By  rotating  the  reference  axis  (figure  3-4)  so 
that  the  magnetic  moment  no  longer  lies  perpendicular  to 
the  field  gradient  in  the  direction  of  force  it  can  be 
shown  a  force  which  has  a  component  in  the  axial  direction 
then  acts  on  the  magnet.  The  cylindrical  symmetry  of  the 
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Figure  3-3 


t 

i 

t 


indicates  current  flow 
indicates  magnetic  field 


'r 


Qualitative  representation  of  the  magnetic 
field  due  to  two  coils  connected  in  a  series 
opposing  configuration. 


Qualitative  representation  of  the  magnetic 
field  due  to  the  two  coils  with  reference 
axis  rotated. 


Figure  3-4 
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suspended  magnet  and  the  calibration  coils  results  in 
the  radial  components  of  force  cancelling.  The  location 
of  the  magnetic  moment  away  from  the  midpoint  of  the 
coil  provides  a  net-force  in  the  axial  direction. 

The  series  opposing  configuration  can  be  posi¬ 
tioned  so  that  the  resulting  magnetic  field  minimizes 
the  e.m.f.  induced  in  the  seismometer  coil.  The  coil 
assembly  with  the  two  coils  at  minimum  separation 
(2.6  cm  centre  to  centre)  was  moved  along  the  seismometer 
until  the  magnitude  of  the  induced  e.m.f.  was  minimized. 
This  was  measured  by  clamping  the  seismometer  and  mea¬ 
suring  the  size  of  the  induced  transient  with  an  oscillo¬ 
scope.  The  position  of  minimum  spike  height  was  found 
to  be  five  and  a  half  centimeters  from  the  seismometer 
base.  The  calibration  coils  are  now  centred  about  the 
fixed  coil  within  the  seismometer. 

The  separation  of  the  two  calibration  coils  was 
then  increased  about  this  point  and  the  displacement  of 
the  magnet  from  its  rest  position  was  measured.  The 
maximum  displacement  (maximum  force  acting  on  magnet) 
was  found  at  a  separation  of  7.6  centimeters.  The  bottom 
of  the  magnet  was  now  entirely  within  the  volume  of  the 
coils.  The  deflection  then  remained  constant  for 
increasing  coil  separation.  The  separation  used  for 
calibration  was  9.6  centimeters.  The  induced  e.m.f.  was 
then  measured  at  2.5  millimeter  intervals  around  the 


s 
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earlier  determined  minimum  with  this  separation  of  the 
coils  (figure  3-5a) .  The  coil  location  is  the  distance 
from  the  bottom  of  the  seismometer  case  to  the  midpoint 
of  the  calibration  coils.  A  new  minimum  was  then  loca¬ 
ted  6  centimeters  above  the  seismometer  base.  This 
change  in  position  was  due  to  a  magnetic  permeable  mass 
(the  magnet)  being  located  between  the  top  coil  and  the 
seismometer  coil.  The  unclamped  response  of  the  seismo¬ 
meter  was  also  measured  and  figure  3-5b  shows  the  varia¬ 
tion  in  magnet  deflection  versus  coil  position. 

3 . 3  Linear  Variable  Differential  Transformer 

An  absolute  calibration  requires  a  knowledge  of 
the  force  which  acts  on  the  suspended  mass  (magnet) .  Thi 
can  be  found  by  measuring  the  axial  displacement,  A&,  of 
the  suspended  mass.  The  applied  force  is  then  kA£  where 
k  is  the  spring  constant  of  the  seismometer.  The  value 
for  is  obtained  using  a  linear  variable  differential 
transformer  (LVDT) .  The  LVDT  used  was  a  500  DC-D  trans¬ 
former  manufactured  by  Schaevitz  engineering.  It  has  a 
sensitivity  of  7.70 ±  0.003  volts  per  centimeter  (figure 
3-6)  and  a  range  of  linearity  of  ±1.2  centimeters. 

The  core  of  the  LVDT  and  its  support  rod  are 
mounted  in  place  of  the  positioning  rod  of  the  seismo¬ 
meter  and  the  transducer  is  mounted  on  top  of  the 
seismometer  over  the  support  rod  (figure  3-7) .  The 
support  rod  is  carefully  and  rigidly  mounted  so  that 
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COIL  LOCATION- cm. 

Figure  3-5a  E.m.f.  induced  in  seismometer  coil  for  a 

one  ampere  current  flowing  through  the 
calibrating  coils  measured  at  various  coil 
locations  (see  text) . 


COIL  LOCATION  -  cm. 

Deflection  of  the  seismometer  mass  in 
millimeters  for  a  one  ampere  current  put 
through  the  calibration  coils  measured  at 
various  coil  locations  (see  text) . 


Figure  3-5b 
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Figure  3-7  The  LVDT  mounted  on  a  Willmore  Mk  II 

seismometer.  The  position  of  calibration 
coils  with  respect  to  the  seismometer  mass, 
clamped  and  unclamped,  is  also  illustrated. 


it  moves  freely  within  the  seismometer  and  LVDT  aperture 
and  therefore  does  not  add  any  damping  to  the  system. 

No  change  is  noted  in  the  amplitude  of  the  output  of 
the  seismometer  when  the  LVDT  is  operating  or  not 
energized.  The  difference  in  the  mass  of  the  suspended 
system  due  to  the  use  of  the  support  and  core  instead 
of  the  positioning  rod  is  24.5  grams,  a  negligible  amount 
when  compared  to  the  4.75  kilograms  of  the  suspended  mass 

The  relationship  for  the  force  on  the  mass  with 
respect  to  current  applied  to  the  coils  can  be  expressed 
as  a  constant.  This  constant  can  be  determined  by 
measuring  the  force  on  the  mass  by  measuring  the  displace 
ment  using  an  LVDT  or  by  originally  calibrating  the 
seismometer  using  one  of  the  methods  in  Chapter  1  and 
using  the  seismometer  to  determine  the  constant.  The 
linearity  of  the  force-current  relation  is  illustrated 
in  figure  3-8.  Force  is  a  constant  times  the  displace¬ 
ment  of  the  seismometer  mass. 

3 . 4  Amplifiers 

To  accurately  record  the  data  the  analog  output 

from  both  the  seismometer  and  LVDT  are  amplified  to 

have  a  dynamic  range  of  ±10  volts.  Both  signals  were 

amplified  using  184  L  operational  amplifiers.  The 

positive  input  was  used  so  that  the  value  of  the  load 

resistor,  RT ,  in  each  case  would  not  affect  the 
L 
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Figure  3-8  Displacement  of  seismometer  mass  for 

various  currents  put  through  the  cali¬ 
bration  coils.  The  position  of  the 
calibration  coils  is  kept  fixed  at  6  cm 
above  the  seismometer  base. 


frequency  response  of  the  amplifier.  In  the  case  of 
the  seismometer  output  it  is  necessary  to  know  the 
frequency  response  of  the  amplifier  used  (figure  3-9a) . 
This  can  be  determined  by  measuring  the  amplitude  and 
phase  response  at  several  frequencies  or  calibrating 
these  values  for  the  amplifier.  The  closed  loop  ampli¬ 
fier  response  is 

z  =  1.0  +  (Rl/R2*x (1.0 + jwRlC) 

where 

Rl  =  4.3  x 10^  ohms 

R2  =  10^  ohms 

C  =  3.36  x  10  9  farads  . 

The  amplitude  response  is  then 

Amp  =  ((Real  z)2  +  (Imag  z)2)  2 

Phz  =  Arctan  (Imag  z/Real  z) 

The  calculated  amplitude  and  phase  responses  are  shown 
in  figure  3-10  along  with  the  measured  values.  The 
close  agreement  of  these  values  show  that  the  calculated 
values  can  be  used  to  transform  the  recorded  data  to  the 
original  frequency  response  of  the  seismometer. 

Only  the  steady  state  value  of  displacement  needs 
to  be  known  so  it  is  only  necessary  to  know  the  D.C. 
gain  of  the  LVDT  amplifier  (figure  3-9b).  The  calculated 
and  measured  values  agree  at  110±1  for  the  D.C.  gain 
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3,300  pf 


Figure  3-9a  The  amplifier  used  to  amplify  the 

output  signal  of  the  seismometer. 


D.C. 

Offset 


+ 15  volts 


Figure  3-9b  The  amplifier  used  to  amplify  the 

output  signal  of  the  LVDT. 
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Figure  3-10  Calculated  phase  and  gain  of  amplifier 

used  on  seismometer  output.  Points 
marked  on  curves  represent  measured 
values . 
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of  the  amplifier.  This  high  gain  as  well  as  the  high 
sensitivity  of  the  LVDT  requires  a  D.C.  offset  on  the 
amplifier  to  prevent  the  amplifier  saturating  when 
the  core  oscillates  about  other  than  the  zero  location 
of  the  LVDT. 

3 . 5  Data  Recording  Sub-system 

The  analog  output  signals  from  the  amplifiers 
were  digitized  so  that  the  data  can  be  manipulated  by 
a  computer.  These  digitized  data  are  then  recorded 
on  a  nine-track,  800  bpi  IBM  compatible  magnetic  tape 
transport.  The  system  of  digitizing  and  recording 
(Allsop  et  al.,  1972)  consists  of  a  modified  model  #120 
Data  Acquisition  System  manufactured  by  Datum  Inc.  and 
a  model  7830  synchronous  tape  transport  made  by  Peri¬ 
pheral  Equipment  Corporations.  In  as  much  as  it  is  a 
synchronous  tape  transport  the  maximum  sampling  rate  is 
governed  by  the  transport  speed.  Since  the  recording 
was  done  on  ten  channels  each  channel  had  a  digitizing 
interval  of  10/2500  =  0.004  seconds.  Each  channel  had 
a  record  length  of  4096  words  and  was  therefore  0.004x4096= 
16.38  seconds  long.  The  response  to  the  step  input  on  each 
record  then  had  adequate  time  to  damp  out  before  the  end 
of  the  record.  The  A  to  D  converter  has  a  resolution  of 
14  bits,  including  sign,  on  a  full  scale  input  of  ±10  volts. 
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Channel  one  was  used  to  record  the  input  step  in  voltage 
measured  across  the  calibration  coils  and  channel  three 
used  to  record  the  output  from  the  LVDT.  The  seismo¬ 
meter  output  is  recorded  on  channels  two,  four,  six, 
eight  and  ten  and  the  remaining  channels  are  grounded. 
The  multiplexer  sampling  interval  is  0.0004  seconds  to 
satisfy  the  synchronous  tape  input  requirements.  If  the 
seismometer  data  are  recorded  on  every  n-th  channel  and 
then  intermixing  the  data  the  sampling  interval  is  n 
times  0.0004.  Recording  on  one  of  the  ten  channels  then 
gives  a  sampling  interval  of  0.004  seconds  and  a  Nyquist 
frequency  (f  =  l/2At  where  At  is  the  sampling  interval) 
of  125  Hz.  Recording  on  every  second  channel  gives  a 
sampling  interval  of  0.0008  seconds  and  a  Nyquist  fre¬ 
quency  of  625  Hz. 
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CHAPTER  4 


THE  PROCESSING  AND  EVALUATION  OF  DATA 
4 . 1  Introduction 

The  purpose  of  calibrating  the  seismometer  is 
to  arrive  at  the  frequency  response  of  the  seismometer- 
log  amplitude  and  phase.  The  frequency  response  of 
the  seismometer  depends  on  the  natural  frequency  (f  ) 
of  the  seismometer,  the  damping  (C)  acting  on  the 
seismometer,  and  the  motor  factor  (A)  of  the  seismometer 
magnet-coil  configuration.  The  frequency  response  is  ob¬ 
tained  by  transforming  the  data  measured  in  the  time  domain 

to  the  frequency  domain,  and  correcting  for  the  frequency 
response  of  the  amplifiers  and  data  recording  subsystem. 


4 . 2  Calculated  Response 


The  expression  for  the  velocity  sensitivity 
(sensitivity  of  the  seismometer  to  an  impulse  in  ground 
velocity)  as  derived  in  Chapter  2  (equation  2.5)  is 

,2 


E  (S)  = 


-AS 


1  +  2CS  +  S' 


=  _  jf 

f 

n 


C  =  damping  factor 
f  =  natural  frequency  in  hertz 
A  =  motor  factor  of  the  seismometer  in  volts 
per  centimeter  per  second. 
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Two  of  the  above  characteristics  are  easily 
altered  on  the  Willmore  Mk  II  seismometer:  the  damp¬ 
ing  factor  and  the  natural  frequency.  Figure  4-1 
shows  the  log  amplitude  and  phase  response  for  E  (S) 
at  several  values  of  ^  -  ranging  from  £=0.3  (the  outer 
curve)  to  £  =  1.0  (the  inner  curve).  The  value  for  A 
is  taken  to  be  1  volt/ (cm/sec)  in  all  cases  and  f  is 
taken  to  be  0.977  Hz,  the  value  found  for  the  seismo¬ 
meter  calibrated.  The  curves  of  the  log  amplitude 

response  have  two  asymptotes  -  one  which  behaves  as 
2 

the  s  term  at  low  frequencies  so  rises  at  40  db  per 
decade  and  one  which  behaves  as  1  at  high  frequencies. 

The  two  asymptotes  intersect  at  the  natural  frequency 

<fn). 

For  £  =  1  (critical  damping)  the  motion  of  the 

suspended  mass  in  the  time  domain  due  to  an  impulse  in 

-a)  t  -uj  t 

ground  velocity  behaves  as  D (e  -  w  t  e  )  where  D 

is  the  magnitude  of  the  inpulse  andu^is  the  natural  fre¬ 
quency  of  the  seismometer  in  radians/sec.  There  is  no 
overshoot  in  the  time  domain.  The  amplitude  response  in 
frequency  domain  is  down  6  db  from  the  horizontal  asymp¬ 
tote  at  the  natural  frequency.  A  more  desirable  response 
is  0.707  critical  damping.  It  is  maximally  flat  at 
higher  frequencies  and  never  exceeds  one  (flat  frequency 
response) .  The  response  is  down  3  db  at  the  natural 
frequency.  A  damping  factor  of  0.5  critical  has  a 


‘ 

-  f 

V 


45 


o 

Q) 

IS) 


Figure  4-1  Calculated  velocity  sensitivity  and  velocity 
phase  of  the  seismometer.  A=  1  volt/ (cm/sec) , 
f n =  .98  hz  and  £  =  0.3,  0.6,  0.707,  0.8  and  1. 
Asymptotes  are  shown  for  the  velocity  res¬ 
ponse  . 
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response  of  1  at  the  natural  frequency  but  exceeds 
the  flat  frequency  response  by  1.25  db  at  1.38  Hz. 
Decreasing  the  damping  made  the  phase  go  through  180° 
more  gradually.  At  low  frequencies  the  seismometer 
output  leads  the  ground  velocity  by  180°  and  at  high 
frequencies  it  is  in  phase  with  the  ground  velocity. 

The  output  always  leads  ground  velocity  by  90°  at  the 
normal  frequency. 

The  damping  for  the  Willmore  Mk  II  seismometer 
is  varied  by  changing  the  load  resistance  as  seen  by 
the  output  of  the  seismometer.  The  value  of  resistance 
needed  to  arrive  at  a  desired  damping  can  be  derived 
from  the  expression 

t 

n  +  1  =  1.7  jj  (4.1) 

where  T  is  the  natural  period  of  the  seismometer  and 

D  is  the  damping  factor  desired.  The  load  resistance 

R  ,  is  related  to  the  resistance  of  the  seismometer 

coil,  R  ,  by  the  relation: 
c 

RL  =  Rc  x  (n)  . 

+ 

Equation  (4.1)  is  derived  in  the  Instruction  Manual 
for  Willmore  Seismometry  Mk  II  put  out  by  Hilger  and 
Watts  Ltd.,  London,  England,  1964.  The  relation 
=  1.6  is  derived  in  Willmore  et  al  (1963). 


n+1 
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For  0.707  critical  damping  at  a  natural  period 
of  one  second  and  a  coil  resistance  of  3.3  kilohms  a 
load  resistance  of  4.29  kilohms  is  needed.  Assumptions 
made  in  the  derivation  of  equation  (4.1)  and  the  pre¬ 
sence  of  additional  resistance  in  the  form  of  the 
seismometer  cable  make  this  value  only  an  approxima¬ 
tion.  The  calibrations  were  done  using  load  resistances 
of  4.11  kilohms,  4.6  kilohms  and  10  kilohms. 

Figure  4-2  shows  the  asymptotes  of  the  log 
amplitude  response  and  the  phase  response  for  inputs 
of  a  step  in  ground  acceleration,  an  impulse  in  ground 
acceleration  (step  in  ground  velocity) ,  and  an  impulse 
in  ground  velocity  (step  in  ground  displacement) .  The 
expression  for  the  seismometer  response  for  a  step  in 
ground  acceleration  is 


E 


step 


(S)  = 


1  +  2<;S  +  s' 


(equation  2.7) 


where  C  is  a  constant.  At  low  frequencies  it  behaves 

2 

as  1  and  at  high  frequencies  it  behaves  as  1/S  .  This 
yields  two  asymptotes  -  one  which  decreases  at  40  db 
per  decade  and  one  which  has  the  value  amplitude  =  1. 
The  two  asymptotes  intersect  at  the  natural  frequency. 
Each  differentiation  of  the  input  signal  causes  a 
rotation  of  +20  db/per  decade  (multiplying  the  ampli¬ 
tude  response  by  the  normalized  frequency  f/fR) . 
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-  response  to  step  in  ground  acceleration 

- response  to  impulse  in  ground  acceleration 

- —response  to  impulse  in  ground  velocity 


Figure  4-2  The  asymptotes  of  the  amplitude  response  and  the 
phase  response  of  the  seismometer  to  various 
inputs  in  ground  motion.  Units  for  the  amplitude 
response  are  volts/ (cm/secj )  for  a  step  in  accel¬ 
eration,  volts/ (cm/sec^ )  for  an  impulse  in  accel¬ 
eration  and  volts/ (cm/sec)  for  an  impulse  in 
ground  velocity. 
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The  phase  response  is  shifted  forward  by  90° 
with  each  differentiation  of  the  input  signal. 

Consider  the  displacement  of  the  suspended  mass  (£) 
and  the  displacement  of  the  seismometer  mass (x) .  For 
low  frequencies  they  are  in  phase  and  for  high  fre¬ 
quencies  £  lags  behind  x  by  180°.  If  we  consider 
the  relative  motion  of  the  suspended  mass  with  respect 
to  the  seismometer  case  (y  =  x  =  £ ,  positive  y  implies 
downward  motion  of  the  suspended  mass  with  respect  to 
the  seismometer  case) ,  the  equation  of  motion  of  the 
suspended  mass  with  no  damping  (equation  2.1)  becomes 

My  +  Ky  =  Mx  . 

This  can  be  rewritten  as 

Ky  =  M*t  . 

The  motion  of  the  suspended  mass  with  respect  to  the 

seismometer  case  (y)  then  is  in  phase  with  the  accel- 

•  • 

eration  of  the  suspended  mass  (£)  which  is  in  phase 
with  the  acceleration  of  the  seismometer  case  (ground 
acceleration)  at  low  frequencies  and  lagging  behind  it 
by  180°  at  high  frequencies.  Since  the  seismometer  is 
a  velocity  transducer  the  output  of  the  seismometer 
(velocity  of  the  suspended  mass  with  respect  to  the 
seismometer  case)  would  then  lead  the  relative  displace¬ 
ment  by  90°.  The  output  would  then  lead  the  ground 
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acceleration  by  90°  at  low  frequencies  and  follow  it 
by  90°  at  high  frequencies  for  an  impulse  in  acceler¬ 
ation.  For  an  impulse  in  velocity  the  seismometer 
output  then  leads  the  ground  velocity  by  180°  at  low 
frequencies  and  is  in  phase  with  ground  velocity  at 
high  frequencies.  The  form  of  input  is  determined  by 
the  signal  put  through  the  calibration  coils.  A  step 
in  current  through  the  coils  is  equivalent  to  a  step 
in  ground  acceleration. 

Varying  the  natural  frequency  (f  )  of  the 
seismometer  changes  the  output  by  moving  it  to  the 
left  or  right  along  the  frequency  scale  as  the  asymp¬ 
totes  always  cross  at  the  natural  frequency.  The  load 
resistance  used  for  a  given  damping  is  a  function  of 
the  natural  frequency  of  the  seismometer  so  a  varia¬ 
tion  in  frequency  without  changing  the  load  resistance 
would  change  the  damping  of  the  seismometer.  The 
natural  frequency  of  the  seismometer  is  varied  by 
varying  the  tension  on  a  "sixth  spoke"  attached  to  the 
magnet  (Willmore  et  al ,  1963)  and  thus  varying  the 
spring  constant,  k,  of  the  seismometer.  The  tension 
is  varied  by  a  screw  adjustment  on  top  of  the  seis¬ 


mometer  case. 


« 
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4 . 3  Mathematical  Treatment  of  the  Data 

The  input  signal,  seismometer  output  and 
LVDT  output  are  recorded  for  four  cases:  no  external 
damping  resistor  (no  electrical  damping)  and  damping 
resistances  of  10  kilohms,  4.6  kilohms  and  4.11  kilohms. 
The  recorded  data  have  a  dynamic  range  of  ±10  volts  in 
14  bits  plus  sign.  A  further  factor  of  4  is  intro¬ 
duced  by  the  program  for  changing  the  data  from  sequen¬ 
tial  storage  into  one  block  per  channel  storage.  Ten 
volts  is  then  represented  by  32,768. 

Any  D.C.  offset  introduced  into  the  seismometer 
output  by  the  amplifier  or  data  recording  subsystem  is 
removed  by  averaging  a  number  of  points  at  the  beginning 
of  the  record  before  the  input  signal  occurs  and  at  the 
end  of  the  record  after  the  signal  has  been  damped  out. 
The  average  is  then  removed  from  the  digitized  data. 

This  does  not  have  to  be  done  with  the  input  signal 
as  it  is  used  only  to  determine  the  initial  condition 
t  (the  starting  point  of  the  response  to  the  input  step) 
or  with  the  LVDT  output  as  this  is  only  used  to  find 
the  steady  state  offset  (Ail)  of  the  suspended  mass. 

Not  removing  the  offset  from  the  seismometer  output 
is  equivalent  to  adding  a  step  in  the  time  domain  over 
the  length  of  the  record  and  would  therefore  offset  the 
transform  of  the  data. 
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The  output  signal  from  the  unterminated  seis¬ 
mometer  (no  load  resistance)  is  used  to  find  the 
natural  frequency  of  the  seismometer.  The  expression 
of  the  output  signal  in  the  time  domain  due  to  a  step 
in  acceleration  is 

AA£u)  -Cw  t  a 

e(t)  =  - 7r~i  e  sin  w  (1-g  )^  t 

(i -rf2  n 

(the  inverse  Laplace  transform  of  equation  2.7). 
Measuring  the  location  and  height  of  the  peaks  yields 

values  for  the  frequency  (w  =  (l-CZ)2u)  )  and  decay 
-Ca^t  °  n 

rate  (e  )  of  the  output  signal.  It  is  then  possi¬ 

ble  to  calculate  the  natural  frequency  (f  =  27Tcon)  and 
the  damping  (C)  of  the  seismometer. 

The  output  signals  from  the  seismometer  with  the 
load  resistors  are  transformed  from  the  time  domain  to 
the  frequency  domain  by  using  z-transforms  (subroutines 
POLYEV,  POLAR  and  DRUM;  Robinson,  1967;  z-transforms 
Kanasewich,  1973).  The  z-transform  operates  in  the  time 
domain  and  can  be  used  at  specific  frequencies.  We 
are  then  able  to  choose  frequencies  at  equal  intervals 
on  a  log  scale,  a  desirable  attribute  as  the  log  ampli¬ 
tude  and  phase  are  both  plotted  vs.  log  frequency. 
Frequencies  are  chosen  so  there  are  twenty  points  per 
decade  between  specified  low  and  high  frequencies. 
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TIME  -  4  milliseconds 

Figure  4-3  The  input  step  digitized  at  4  millisecond 
intervals. 


Figure  4-4a  The  phase  response  of  the  seismometer  with 

a  10  kilohm  load  resistor  digitized  at 
intervals  of  4  millisecond  and  having  a 
Nyquist  frequency  of  125  Hz.  The  starting 
point  of  transformation,  t  ,  is  taken  to  be 

tf. 
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Figure  4-4  b-d  The  phase  response  of  the  seismometer  with 

a  Nyquist  frequency  of  125  hz  using  values 
of  t  equal  to  t^ ,  t^  and  t^ . 
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The  starting  point  of  the  transformation,  t 
is  determined  for  the  input  signal.  Figure  4-3  shows 
the  recorded  input  step.  The  dependence  of  phase  on 
the  point  chosen  for  t  is  shown  in  figure  4-4.  The 
phase  response  varies  by  180°  at  the  Nyquist  frequency 
(125  Hz)  for  every  successive  point  chosen  as  t  . 

Figure  4.4c  shows  the  phase  response  which  uses  the 
first  point  after  the  step  occurred  as  t  ;  the  value 
used  in  the  calibration. 

The  amplitude  of  the  transformed  data  at  each 
frequency  was  divided  by  the  calculated  gain  of  the 
seismometer  amplifier  at  that  frequency  and  by  the 
calculated  value  for  A£  (determined  by  dividing  the 
steady  state  change  in  volts  of  the  LVDT  output;  by  the 
gain  of  the  LVDT  amplifier,  110  ±  1;  and  by  the  sensiti¬ 
vity  of  the  LVDT;  7.70  ±0.003  volts/cm).  The  amplitude 

3 

is  now  expressed  in  volts/ (cm/sec  ).  This  is  the  seis¬ 
mometer  response  to  a  step  in  ground  acceleration. 

The  amplitude  is  then  rotated  by  multiplying  each  point 

of  the  frequency  response  by  the  normalized  frequency 

2 

response  at  that  point  squared  (  (f/f  )  )  to  arrive  at 

the  velocity  sensitivity. 

The  phase  response  is  corrected  for  the  calcu¬ 
lated  phase  change  due  to  the  seismometer  amplifier 
and  delay  introduced  since  all  channels  are  not  digi¬ 
tized  simultaneously.  This  introduces  a  continuous 
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phase  shift  which  reaches  360°  at  2500  Hz  (2500  = 

1/.0004  where  .0004  is  the  delay  between  channels) 
or  a  shift  of  18°  at  125  Hz.  The  phase  shift  is 
changed  to  one  due  to  an  impulse  in  velocity  by 
shifting  it  forward  by  180°.  The  log  amplitude 
and  phase  are  then  plotted  against  log  frequency. 

4 . 4  Observed  Results 

Calibration  runs  were  done  on  a  Willmore  Mk  II 
seismometer,  serial  number  56450/229518.  The  runs  are 
done  using  step  and  sine  wave  inputs  in  current  through 
the  calibrating  coils,  and  using  the  Maxwell  bridge. 

The  data  from  the  step  input  are  digitized  and  pro¬ 
cessed  as  in  section  4.3.  The  output  from  the  other 
two  methods  is  read  off  an  oscilloscope. 

a)  A  Step  in  Input  through  the  Calibration  Coils 

A  step  in  current  through  the  calibration  coils 
is  used  to  apply  a  force  to  the  suspended  magnet  within 
the  seismometer.  This  step  in  magnetic  field  is  equi¬ 
valent  to  a  step  in  the  inertial  frame  which  is  equi¬ 
valent  to  a  step  in  ground  acceleration.  Figure  4-5 
shows  the  log  amplitude  response  observed  due  to  a  step 
in  ground  acceleration  (step  acc.  sens.)  for  a  4.1  kilohm 
load  resistor  and  the  velocity  sensitivity  of  the  seismo¬ 
meter  calculated  from  it.  The  digitizing  interval  used 
was  0.004  seconds  yielding  a  Nyquist  frequency  of  125  Hz 


, 
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Figure  4-5  The  seismometer  response  to  a  step  in  ground 
acceleration  and  the  calculated  velocity  sen¬ 
sitivity  for  a  seismometer  with  a  load  resis¬ 
tance  of  4.1  kilohms.  The  data  is  digitized  at 
an  interval  of  4  milliseconds.  The  horizontal 
asymptote  has  a  value  of  2.9. 
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The  asymptotes  behave  as  expected  in  each  case  -  one 
asymptote  of  the  step  acc.  sens,  plot  decreases  at 
40  db/decade  at  the  higher  frequencies  while  at  lower 
frequencies  the  asymptote  is  level  at  about  2.9 
volts/ (cm/sec  ).  The  velocity  sensitivity  plot  has  an 
asymptote  with  a  slope  of  40  db/decade  at  lower  fre¬ 
quencies  and  a  horizontal  asymptote  at  2.9  volt/ (cm/sec) 
at  higher  frequencies.  Both  asymptotes  in  each  case 
intersect  at  the  natural  frequency  of  the  seismometer 
(0.98  Hz) . 

An  increase  in  the  observed  values  from  the  ex¬ 
pected  values  at  high  frequencies  can  be  seen  on  each  of 
the  two  graphs  in  figure  4-5.  As  the  corner  of  the 
seismometer  amplifier  was  at  100  Hz  and  the  Nyquist 
frequency  at  which  the  data  was  digitized  was  125  Hz 
it  was  believed  that  the  observed  increase  may  be  due 
to  alliasing.  To  lessen  this  effect  the  data  was  digi¬ 
tized  at  an  interval  of  0.0008  seconds  giving  a  Nyquist 
frequency  of  625  Hz.  Figures  4-6  to  4-8  show  the  log 
amplitude  curves  for  the  seismometer  output  signal  di¬ 
gitized  at  this  interval.  Comparing  the  step  acceleration 
sensitivity  and  velocity  sensitivity  for  a  seismometer  with 
a  damping  resistor  of  4.1  kilohms  (figure  4-8)  with  the 
step  acceleration  sensitivity  and  velocity  sensitivity  of 
the  seismometer  with  the  same  damping  but  digitized  at 
0.004  second  interval  (figure  4-5), shows  no  significant  change 
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FREQUENCY  -  hertz 

Figure  4-6  The  seismometer  response  to  a  step  in  ground 
acceleration  and  the  calculated  velocity  senr 
sitivity  of  a  seismometer  with  a  damping  re¬ 
sistance  of  10  K.  The  digitizing  interval.  At, 
is  .0008  seconds.  The  asymptotes  are  indicated 
by  a  dashed  line. 


VEL.  SENS.  —  voltsAcm/sec)  SfEP  ACC.  SENS.  -  volts /{cm/sec3) 
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Figure  4-7  The  seismometer  response  to  a  step  in  ground 
acceleration  and  the  calculated  velocity  sen¬ 
sitivity  of  the  seismometer,  At  =  .0008 
seconds . 
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Fiaure  4-8  The  seismometer  response  to  a  step  in  ground 
acceleration  and  the  calculated  velocity 
sensitivity  for  a  seismometer  with  a  4.1  K 
damping  resistor  At  =  .0008. 
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in  the  outputs.  It  must  then  be  concluded  that  the 
observed  increase  is  not  primarily  due  to  alliasing 
but  rather  to  some  other  source. 

A  possible  source  can  be  seen  by  examining 

the  size  of  this  increase.  The  round  off  of  the  data 

in  digitizing  gives  us  a  possible  error  of  the  size  of 

half  a  bit  on  the  A-D  converter.  This  would  result  in 

-3  3 

an  increase  of  1.98  x 10  volts/ (cm/sec  )  in  the  step 
acceleration  sensitivity  values.  When  the  values  are 
differentiated  to  obtain  the  velocity  sensitivity  this 
introduces  an  error  of  .8  volts/ (cm/sec)  at  20  Hz,  3.2 
volts/ (cm/sec)  at  40  Hz  and  12.8  volts/ (cm/sec)  at 
80  Hz.  Values  less  than  this  are  all  that  are  needed 
to  account  for  the  noted  increase  in  the  velocity  sensi¬ 
tivity.  The  observed  step  acceleration  sensitivity  shows 

-3  3 

an  increase  of  about  1.4  x  10  volts/ (cm/sec  )  over  the 
expected  response.  This  falls  well  within  the  quantiza¬ 
tion  error  of  the  A-D  converter  (Plesinger,  1971) . 

The  velocity  sensitivity  of  the  seismometer 
with  a  damping  resistor  of  10  kilhms  shows  a 
horizontal  asymptote  of  4.05  volts/ (cm/sec)  and  a 
damping  of  0.35  critical  (calculated  from  the  graph 
by  reading  off  the  response  of  the  seismometer  at  the 
natural  frequency) .  The  motor  factor,  A,  of  the  seis¬ 
mometer  is  not  the  value  of  the  horizontal  asymptote 
as  the  output  signal  was  attenuated  by  the  load 
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resistance  R 
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V 


V 


obs  RT  +  R  'out 
L  s 


where  is  the  observed  output,  VQUt  is  the  output 

from  the  seismometer,  and  R  is  the  resistance  of  the 

s 

coil  within  the  seismometer.  The  motor  factor  of  the 
seismometer  is  then  R  +RT/RT  times  the  value  of  the 
horizontal  asymptote.  The  measured  seismometer  coil 
resistance  is  3.68  kilohms.  This  yields  a  motor  factor 
of  5.53  volts/ (cm/sec)  for  the  10  kilohm  load  resistor. 
The  4.6  kilohm  load  resistor  shows  a  damping  of  0.45 
critical  and  horizontal  asymptote  of  3.05  volts/ (cm/sec) 
yielding  a  motor  factor  of  5.49  volts/ (cm/sec) .  A 
damping  of  0.62  critical  and  horizontal  asymptote  of 
2.9  volts/ (cm/sec)  are  observed  for  the  4.1  kilohm  load 
resistor  giving  a  motor  factor  of  5.52  volts/ (cm/sec) ; 
this  yields  an  average  motor  factor  of  5.51 ±  0.01  volts/ 
(cm/sec)  for  the  seismometer.  The  observed  damping 
indicated  that  0.707  critical  damping  requires  a  load 
resistance  slightly  less  than  4.1  kilohms  and  shows  the 
error  in  the  method  of  calculating  the  damping  resistance 
Figure  4-9  shows  the  velocity  phase  response  of 
the  seismometer  with  different  damping  resistance.  The 
phase  behaves  as  expected,  leading  the  ground  velocity 
by  180°  at  low  frequencies  and  being  in  phase  with  it 
at  high  frequencies.  The  phase  response  above  1.8  Hz 
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Figure  4-9  Velocity  phase  response  of  the  seismo¬ 
meter  for  various  load  resistances. 


65 


is  not  valid  because  of  the  decreased  signal  to  noise 
ratio . 

b)  A  Sinewave  Input  through  the  Calibrating  Coils 

A  calibration  run  was  done  by  putting  a  sine- 
wave  through  the  coil  at  given  frequencies  and  measur¬ 
ing  the  frequency  response  by  measuring  the  height  of 
the  sinewave  observed  in  the  output.  The  phase  was 
measured  using  a  variable  phase  signal  generator. 
Comparison  of  the  two  frequency  responses  would  show 
any  error  in  the  step  calibration  method  due  to  the 
assumed  shape  of  the  step,  the  digitization  of  the  data, 
the  transformation  of  the  data  or  alliasing.  Reading 
the  peak  to  peak  values  of  the  output  sinewave  and 
dividing  by  the  steady  state  deflection  of  the  magnet 
yield  the  seismometer  sensitivity  to  an  impulse  in 
ground  acceleration  (acc.  sens.).  The  steady  state 
deflection  of  the  magnet  is  obtained  from  a  graph  of 
the  deflection  vs.  frequency  on  a  log  scale  (figure 
4-10).  The  steady  state  deflection  for  the  seismometer 
mass  is  0.0127  cm  for  a  10  kilohm  load  resistance, 

0.0148  cm  with  a  4.6  kilohm  load  resistance  and  0.0145 
cm  with  a  4.1  kilohm  load  resistance.  The  acceleration 
sensitivity  of  the  seismometer  calculated  using  these  va¬ 
lues  is  shown  in  figure  4-11  and  the  velocity  sensitivity 
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Ficrure  4-11  Response  of  the  seismometer  to  an  impulse 

in  ground  acceleration.  Asymptotes  are 
indicated  by  a  dashed  line.  The  natural 
frequency  is  .98  hz. 
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calculated  from  the  acceleration  sensitivity  is  shown 
in  figure  4-12.  The  velocity  sensitivity  shows  hori¬ 
zontal  asymptotes  of  4.10  volts/ (cm/sec)  for  a  load 
resistance  of  10  kilohms,  3.10  volts/ (cm/sec)  for  a 
load  resistance  of  4.6  kilohms,  and  2.9  volts/ (cm/sec) 
for  a  load  resistance  of  4.1  kilohms.  These  values 
agree  with  those  arrived  at  from  the  step  input.  An 
overlay  of  the  two  sets  of  curves  shows  an  exact 
agreement  within  the  useable  frequency  band.  The 
measured  phase  is  shifted  to  represent  the  phase  between 
the  measured  signal  and  ground  velocity  by  adding  90° 
(figure  4-13)  .  This  too  is  in  agreement  with  the  cal¬ 
culated  values. 

c)  Maxwell  Bridge  Technique 

For  comparison  purposes  a  calibration  run  was 
made  on  the  seismometer  using  the  Maxwell  bridge 
technique  (Chapter  1  section  1.2d).  Load  resistances 
of  10  kilohms  and  4.1  kilohms  were  used.  The  results 
were  evaluated  using  a  program  written  by  Dr.  E.R. 
Kanasewich.  The  velocity  sensitivity  arrived  at  by 
this  method  is  shown  in  figure  4-14.  The  4.1  kilohm 
resistor  yielded  a  horizontal  asymptote  of  2.8-volts/ 
(cm/sec)  and  a  damping  of  0.63  critical.  The  motor 
factor  of  the  seismometer  calculated  by  this  technique 
was  5.72  volts/ (cm/sec) .  The  10  kilohm  resistor  had 
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FREQUENCY- hertz 

Figure  4-12  The  velocity  sensitivity  of  the  seismometer 

with  various  load  resistances  for  a  sine-wave 
input.  The  asymptotes  are  indicated  by  dashed 
lines . 
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Figure  4-13  Velocity  phase  response  for  various 

damping  resistances  measured  using  a 
variable  phase  signal  generator  and 
a  sine-wave  input. 


VEL.  SENS. -vo!ts/(cm/sec)  VEL,  SENS  -volts/Ccm/sec) 
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FREQUENCY  -  HERTZ 

Figure  4-14  The  velocity  sensitivity  of  the  seismometer 

arrived  at  using  the  Maxwell  bridge  calibra¬ 
tion  technique  (Willmore,  1959). 
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a  horizontal  asymptote  of  3.4  volts/ (cm/sec) ,  a 
damping  of  0.36  critical  and  a  motor  factor  of  6.01 
volts/ (cm/sec) .  The  damping  change  due  to  the  use 
of  the  Maxwell  bridge  is  that  the  load  resistance 
as  seen  by  the  seismometer  was  not  the  load  resistance 
put  on  the  amplifier  but  rather  the  load  resistance  of 
the  amplifier  in  parallel  with  a  variable  resistance 
of  about  200  kilohms  in  size.  The  load  resistance 
as  seen  by  the  seismometer  is  then  less  than  the  values 
used  (1.0  kilohms  and  4.1  kilohms)  and  so  the  signal  is 
more  heavily  damped  than  the  cases  using  the  calibration 
coils.  The  difference  in  the  values  of  the  motor  factor 
of  the  seismometer  as  measured  by  the  calibration  coils 
(5.51  0.01  volts/ (cm/sec) )  and  the  Maxwell  bridge 

technique  (5.72  volts/ (cm/sec)  and  6.01  volts/ (cm/sec) ) 
may  be  due  to  the  fact  that  the  Maxwell  bridge  calibra¬ 
tion  was  carried  out  over  a  month  after  the  calibration 
with  the  calibrating  coils  and  the  motor  factor  of  the 
seismometer  may  have  changed  in  this  time  but  the 
disagreement  in  the  motor  factor  for  the  two  runs  done 
on  the  Maxwell  bridge  using  two  different  damping  resis¬ 
tances  indicates  a  possible  source  of  error  is  the 
reading  taken  by  this  method  or  by  the  method  itself. 
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4 . 5  Conclusion 

A  step  in  current  through  the  calibrating  coils 
provides  a  rapid,  accurate  method  of  calibrating  any 
moving  magnet  seismometer.  Though  the  calibration  run 
was  done  only  on  a  vertical  seismometer  neither  the 
direction  nor  magnitude  of  the  applied  force  depends  on 
gravity.  The  calibration  technique  is  then  independent 
of  the  orientation  of  the  seismometer  and  calibrating 
coils  with  respect  to  the  earth. 

Alliasing  and  digitizing  errors  involved  in  the 
transformation  of  the  data  into  the  frequency  domain 
can  be  minimized  by  using  a  least  squares  fit  on  the 
data  to  determine  the  natural  frequency,  damping  and 
value  of  the  horizontal  asymptote  of  the  seismometer 
response.  The  velocity  sensitivity  can  be  found  by 
substituting  these  values  into  the  correct  expression. 

Whether  the  velocity  sensitivity  is  found  using 
a  least  squares  fit  or  transforming  the  data  into  the 
frequency  domain  this  method  is  rapid  as  it  only 
requires  the  input  of  one  signal,  a  square  wave  and 
the  digitizing  of  the  seismometer  output.  The  sine- 
wave  technique  requires  readings  at  a  large  number  of 
frequencies  to  ensure  accuracy  of  data  so  is  a  time 
consuming  method  and  therefore  not  easily  applicable 
to  frequent  calibration  of  seismometers.  The  simplicity 
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of  the  step  input  is  ideal  in  that  it  can  be  generated 
by  a  clock  D-A  converter  and  so  no  signal  generator  is 
needed. 

The  linearity  of  the  deflection  of  the  seismo¬ 
meter  magnet  with  respect  to  input  current  makes  it 
unnecessary  to  measure  the  magnet  deflection  once  the 
coils  are  mounted  on  the  seismometer  case  and  an  initial 
calibration  is  done.  A  remeasurement  of  deflection  is 
only  necessary  if  a  sizable  change  is  noticed  in  the 
motor  factor  of  the  seismometer.  Only  one  linear 
variable  differential  transformer  is  needed  to  cali¬ 
brate  all  seismometers. 

As  no  contact  with  the  seismometer  is  needed 
this  technique  is  easily  used  for  frequent  field  cali¬ 
brations  where  the  seismometer  may  be  buried  to  reduce 
the  effects  of  pressure  and  temperature  variation. 

This  method  then  has  advantages  of  the  traditional 
calibration  techniques  such  as  Willmore's  use  of  the 
Maxwell  bridge  and  the  use  of  shake  tables. 
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